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We show that electrons on hquid hehum display intrinsic bistability of resonant inter-subband 
absorption. The bistability occurs for comparatively weak microwave power. The underlying giant 
nonlinearity of the many-electron response results from the interplay of the strong short-range 
electron correlations, the long relaxation time, and the multi-subband character of the electron 
energy spectrum. 
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Electrons on liquid helium provide a unique tool 
for studying correlation effects in two- dimensions (2D). 
The ratio of the characteristic Coulomb energy to the 
in-plane kinetic energy, the plasma parameter V = 
e^iiinsY^'^ /ksT {ug is the electron surface density), can 
vary by orders of magnitude, from F <C 1 where the elec- 
tron system is a weakly-interacting gas to F > 130 where 
it is a Wigner crystal (TJ |2j . In the broad range 1 < F < 
130 the electron system is a correlated 2D electron liquid, 
with unusual and sometimes counter- intuitive properties 
of classical and quantum 2D magneto-transport and 
activated and tunneling escape from the surface [4]. 

The effect of correlations in an electron liquid was also 
seen in weak-field spectroscopy of transitions to excited 
subbands of motion along the helium surface [5 . In a 
sense, it is a counter-part, for a strongly correlated sys- 
tem, of the depolarization effect in inter-subband absorp- 
tion in semiconductor heterostructures with high elec- 
tron densities [6 . For semiconductors, of much interest 
are both correlation transport effects and nonlinear 
optical effects related to radiation-induced population of 
excited subbands [8 . Long sought has been optical bista- 
bility in intersubband absorption [9 . 

For electrons on helium, optical nonlinearity should be 
strong, since the electron relaxation time is unusually 
long, reaching ~ 10~^ s for T ^ 0.1 K, and saturation of 
resonant inter-subband absorption has been indeed seen 
[To] . The interplay of the long relaxation time and strong 
spatial correlations provides a qualitatively new nonlin- 
earity mechanism and should lead to new resonant ef- 
fects. Such effects are indeed found in this paper. 

Our central result is the first, to the best of our knowl- 
edge, direct experimental observation and a theory of 
the bistability of resonant inter-subband absorption in 
a correlated electron system. The bistability is due to 
the correlation-induced strong dependence of the inter- 
subband transition frequencies of an electron on the state 
of other electrons. An additional interest in this depen- 
dence comes from the proposals of quantum computing 



with electrons on helium [TT] , as it provides a mechanism 
of two-qubit gate operations. 

We study electrons on liquid ^He. The experimental 
setup is similar to that previously described [12 . Elec- 
trons are confined on a helium surface between two circu- 
lar parallel electrodes in an asymmetric potential formed 
by the barrier on the surface, the image force, and an elec- 
tric field Ej_ normal to the surface from the voltage on 
the electrodes [13j. The quantized energies of ID motion 
in this potential ei {I = 1,2...) give the inter-subband 
transition frequencies uom = (e^/ — ei)/fi in the single- 
electron approximation. We resonantly excite 1^2 
inter-subband transitions with microwave (MW) radia- 
tion {uo/ 271 = 104.5 GHz) by tuning the frequency UO21 
with the field through the linear Stark shift. 

The resonant MW response is observed by measuring 
the low-frequency (100 kHz) in-plane magnetoconductiv- 
ity (J which depends on the electron temperature Tg. 
We apply a classically strong magnetic field B ~ 350 G 
normal to the surface and use a Corbino disk that con- 
stitutes the top electrode. In these conditions, because 
of the strong electron correlations, cFxx ^ c^UqI' / {muoD^ 
where uOc = eB/mc and v is the single-electron scatter- 
ing rate with no magnetic field [3 . This rate, in turn, 
is determined by Tg [14 . Thus the variation of cjxx re- 
flects the degree of electron heating caused by resonant 
MW absorption and allows us to determine Tg from the 
relative change of the electron scattering rate v [12 . 

The underlying physical picture is based on the hi- 
erarchy of relaxation processes in the resonantly modu- 
lated correlated electron system [15]. The fastest pro- 
cess is the exchange of in-plane momentum between the 
electrons. It occurs over the typical time oo~^ ^ where 

uOp = (27re^ns^^/m)^/^ is the electron plasma frequency. 
The associated exchange of in-plane energy leads to es- 
tablishing an in-plane electron temperature Tg, which is 
the same for all occupied subbands. Next comes the rate 
V of quasi-elastic momentum scattering by helium va- 
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por atoms or by surface waves on helium, ripplons. This 
scattering is short-range, and therefore v characterizes 
not only the in-plane momentum relaxation rate, but 
also the rate of elastic inter-subband transitions. Such 
transitions lead to thermal distribution over subbands 
characterized by the same temperature Tg. The slowest 
process is the electron energy relaxation, which is due to 
inelastic processes, including one- and two-ripplon scat- 
tering or scattering by bulk excitations in helium pT] . 

According to the above picture, when a MW field res- 
onantly excites 1^2 transitions, the energy goes first 
to heating up the electron system. As a result electrons 
populate higher subbands and the lowest subband is de- 
pleted, which leads to the decrease of resonant absorp- 
tion. Such absorption bleaching occurs for much smaller 
radiation intensity than necessary for absorption satura- 
tion in a two-subband model [121 ES] • 

Figure [l] shows the examples of the measured resonant 
response for = 4.0 x 10^ cm~^ at T = 0.4 K, plotted 
as /S.crxx/^xxi where /S.crxx is the change of the electron 
conductivity when MW radiation is applied. Each trace 
was obtained by slowly increasing E±_ to drive electrons 
through resonance, while maintaining the input radiation 
power P at a fixed level. At T = 0.4 K the scattering 
rate v is dominated by collisions with He vapor atoms 
and monotonically increases with Tg [T2j [14]. Corre- 
spondingly, upon sweeping E±_ , cjxx increases and attains 
maximum at the resonance. 

It is seen from Fig. [l] that as the power increases, 
the resonance shifts towards lower field values. This 
means that the transition frequency UO21 increases rela- 
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FIG. 1: Conductivity change Aaxx/cTx 



E± for the surface density ris = 4.0 x 10^ cm ^ 



XX vs. the tuning field 
~ at T = 0.4 K. 

Solid lines are traces obtained in an upward sweep. Corre- 
sponding to their increasing amplitude, the traces were taken 
at MW power P = 40, 130, 380 and 1000 /iW. The dashed 
line is the trace taken in a downward sweep at P = 1000 /iW 
(the sweeping directions are indicated by arrows). 



tive to its value for electrons at Tg ~ T. The frequency 
shift Acoc obtained from the experiment is plotted in 
Fig. [2] as a function of T^. The conversion factor to 
frequency (the slope of a;2i/27r versus E±) was deter- 
mined experimentally by varying the frequency of the 
MW source and recording the corresponding shift of the 
resonance along the E± axis. This factor was found to 
be 0.64 ± 0.03 GHz-cm/V for the resonance in Fig. [l] 
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FIG. 2: Frequency shift Acjc/27r vs. electron temperature 
Te obtained from the measured resonant response (squares). 
The solid line is given by Eq. ^ with P=8.91. 

The major cause of the experimental frequency shift 
is the electron-electron interaction. The corresponding 
part of the Hamiltonian in a correlated system is 



^f/K 



(1) 



4 ^ ^n^ra 

where and are the lateral and normal to the helium 
surface electron coordinates, respectively. We have taken 
into account that \zn — Zm\ <C Ir^ — r^|. The distance 
— is given by the localization length of the quan- 
tized states |/) of motion normal to the surface, which 
is ~ 10 nm [13], whereas the inter-electron distance is 
^ 1 /im for typical Ug and Tg; note that (i^ee^) would 
diverge if there were no spatial electron correlations. 

Assuming that the lateral positions of the electrons and 
their distribution over the subbands are uncorrelated, to 
first order in i^ie^ one obtains for the frequency shift [15] 



ACJ21 = 

-2(zii 



2h 



Z22) 



ziipii ^2\zi2f{pii - P22) , (2) 



where = {l\A\V); the matrix element pu is the 

fractional occupation of subband and F is related 
to the electron pair-correlation function g{r) as F = 



-1/2 



Jr- 



g{r)dY. Equation \2\ differs from the cor- 



responding expression in Ref. [5 even for pu = Si^i, as 
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assumed there. Evidence of a transition frequency shift 
that depends on subband occupations was found by Glas- 
son et al [16]. 

If the many-electron system is described by the Boltz- 
mann distribution with temperature Tg, the transition 
frequency shift, Eq. is a function of Tg. The depen- 
dence of Au2i on Tg comes from the occupations pu and 
the factor F that is determined by the plasma parameter 
r calculated for T = T^. For an electron crystal with 
triangular lattice F ^ 8.91. The experiment [5 showed 
that F does not increase significantly for 10 ^ F < 100. 
This agrees with the Monte-Carlo simulations For 
experimental parameters and T of Fig. [T] we have 
F ~ 50 for Tg ~ T. Therefore, to estimate Aa;2i we used 
F = 8.91. The result for the electron- heating induced 
frequency shift AuJc = Auj2i{Te) — Au2i{T) is shown by 
the solid line in Fig. |2j The observed shift is larger by 
about 30% than the theoretical estimate. 

The electron-electron interaction also affects the line 
shape of the resonance as seen in Fig. [l] As the MW 
power P increases, the resonance line becomes asymmet- 
ric with a steeper low-field side. At high powers the line 
shape changes dramatically, as illustrated in Fig.[l]by ex- 
perimental traces taken at P = 1000 jaW. In an upward 
sweep, the low-field side becomes an abrupt jump that 
corresponds to a sudden increase of Tg. In a downward 
sweep, the response coincides with the upward- sweep 
curve on the high-field side of the resonance. However, it 
takes on substantially higher values on the low-field side. 
This unusual behavior becomes even more pronounced 
at lower temperatures. Figure [3] shows a high power re- 
sponse for Us = 8.0 X 10^ cm-2 and T = 0.2 K. At this 
temperature, the scattering rate u (therefore axx) is de- 
termined by the interaction of electrons with ripplons and 
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FIG. 3: Resonant response for Us = 8.0 X 10^ cm"^ at 
T = 0.2 K obtained in an upward (open squares) and down- 
ward (solid squares) sweeps with P = 870 /iW. The sweeping 
directions are indicated by arrows. 



decreases with Tg \TT . Correspondingly, the response is 
plotted as —Aa^xl^xx- 

The hysteresis of resonant response indicate the on- 
set of bistability, i.e. co-existence of two stable regimes 
characterized by different values of Tg. This effect can be 
explained by considering the balance between the heat- 
ing of electrons by the MW radiation and cooling by the 
thermal bath (liquid helium). The energy balance equa- 
tion can be written as [12l [15] 

hujT (pii - P22) = VE^B (Tg - T) . (3) 

Here, ve is the electron energy relaxation rate and 
r is the rate of MW-induced transitions, r = 
0.5J^^7/ [(^ + Acj2i)^ +7^], where VLr = eE^zx2l^ is 
the Rabi frequency (T^o is the MW electric field), 7 is 
the half- width of the absorption line, and ^ = {jJ2\ — ^ 
is (minus) the MW frequency detuning from the single- 
electron resonance. 

The solution of Eq. (|3| gives the stationary temper- 
ature of the electron system. A way of obtaining this 
solution graphically is illustrated in Fig. [4j The solution 
is defined by the intersection of a solid line and a dashed 
line, which represent the cooling rate and the heating 
rate, respectively. For the numerical evaluation, we as- 
sumed that is determined by collisions with He vapor 
atoms. However, we note that for ^He at T = 0.4 K the 
energy relaxation rate can be affected by the interaction 
of electrons with short-wavelength ripplons, cf. [iT]. 
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FIG. 4: Graphical solutions of the energy balance equation 
([3]) for Us = 4.0 X 10'^ cm"^ and T = 0.4 K. The solid line 
shows the cooling rate VE^uiTe — T). The dashed lines show 
the heating rate Tiuori^pn — P22) calculated for three frequency 
detunings CO21 — a; = ^1, ^bi and 62 (-1.15, -0.85 and -0.7 GHz, 
respectively). The other parameters are Q.r/2ti = 0.01 GHz, 
7/27r = 0.2 GHz and ijj/2n = 104.5 GHz. The circles mark 
two stable solutions of Eq. Q for detuning ^bi- Inset: Te vs. 
5 — U21 — 00 calculated from Eq. ([3| for the same ns,T, Qr, 
7, and 00. The high- and low-Te branches are plotted by the 
solid lines, the dashed arrows indicate inter-branch switching. 

Important for the onset of bistability is that the spec- 
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tral peak of resonant absorption rate r{Lj) is narrow and 
that its position depends on Tg in terms of Acj2i(^e)- 
This leads to multiple solutions of Eq. (|3| in a certain 
frequency range for sufficiently large MW power, as seen 
from Fig. |4] [15 . In particular, there exist two stable 
states (indicated by circles in Fig.|4|: the one with com- 
paratively high Tg and another one with comparatively 
low Tg. The high-Tg state ceases to exist for 6 < di^ 
while the low-Tg state is not possible for S > S2. The 
boundaries 5i and S2 of the bistability interval depend 
on radiation power and shift towards negative values with 
increasing ftR. 

This theoretical model explains well the experimen- 
tally observed behavior. The calculated Tg as a function 
of S is shown in the inset of Fig. |4] Two stable solutions 
are plotted by the solid lines. As we increase S starting 
from a large negative value, i.e. approach the resonance 
by sweeping E± from a low- field side, the system follows 
the low-Tg branch as long as this solution exists. At the 
end of the branch, the system undergoes a transition to 
a new equilibrium state with higher Tg, as indicated by 
the upward arrow. Correspondingly, in the experiment 
we expect to observe an abrupt jump of the resonance re- 
sponse. On the other hand, when we decrease S starting 
from a large positive value, i.e. sweep through the res- 
onance from a high-T^^ side, the system remains on the 
high-Tg branch until it ends. Then the system switches 
to the small- Tg branch, as shown by the downward arrow. 
In the experiment, it could be seen as a sudden decrease 
of the resonant response. 

Our model assumes that the system is spatially uni- 
form and the absorption linewidth 7 is independent of Tg. 
Heating-induced broadening of the absorption line might 
account for the difference in the response for T = 0.4 K 
and 0.2 K (see Figs, [l] and [s]). For T = 0.4 K and for high 
Tg the linewidth 7 becomes comparable to the inhomoge- 
neous broadening, extending the frequency range where 
the system stays on the high-Tg branch. At T = 0.2 K 
this does not happen, 7 is much smaller, and the shape 
of the observed response in closer to the calculated one. 
A more detailed discussion will be provided elsewhere. 

In summary, we have directly observed the long-sought 
intrinsic optical bistability in resonant inter-subband ab- 
sorption of a quasi-2D electron liquid. The effect is due 
to the electron-electron interaction, which leads to strong 
spatial correlations in the electron liquid and as a conse- 
quence, to a strong shift of the inter-subband absorption 
line with varying distribution over the electron states. 
The bistability was studied by measuring the electron 
magnetoconductivity. The proposed theoretical model is 
in good agreement with the experimental results. 

We acknowledge valuable discussions with P. M. Platz- 
man and D. Ryvkine. DK and KK were supported in 
part by the MEXT, Grant-in- Aids for Scientific Research, 
MID was supported in part by the NSF, grant EMT/QIS 
0829854, MJL was supported in part by the EPSRC. 



* E-mail: konstantinov@riken.jp' 
[1] C. C. Grimes and G. Adams, Phys. Rew. Lett. 42 795 

(1979); D. S. Fisher, B. 1. Halperin, and P. M. Platzman, 

Phys. Rev. Lett. 42, 798 (1979). 
[2] Electrons on Helium and Other Cryogenic Substrates, 

edited by E.Y. Andrei (Kluwer Academic, Dordrecht, 

1997); for an introductory review see W. F. Vinen and 

A. J. Dahm, Phys. Today 40, 43 (1987). 
[3] M. I. Dykman and L. S. Khazan, JETP 50, 747 (1979); 

M. I. Dykman, M. J. Lea, P. Fozooni, and J. Frost, Phys. 

Rev. Lett. 70, 3975 (1993); Yu P. Monarkha, S. Ito, 

K. Shirahama, and K. Kono, Phys. Rev. Lett. 78 2445 

(1997) ; M. J. Lea and M. I. Dykman, Physica B 251, 628 

(1998) . 

[4] Y. lye, K. Kono, K. Kojita, and W. Sasaki, J. Low 
Temp. Phys. 38, 293 (1980); G.F. Saville, J.M. Good- 
kind, and P.M. Platzman, Phys. Rev. Lett. 70, 1517 
(1993); L. Menna, S. Yiicel, and E.Y. Andrei, ibid. 70, 
2154 (1993); M. 1. Dykman, T. Sharpee, and P. M. Platz- 
man, ibid. 86, 2408 (2001). 

[5] D.K. Lambert and P.L. Richards, Phys. Rev. Lett. 44, 
1427 (1980). 

[6] S. J. Allen, D. C. Tsui, and B. Vinter, Solid State Com- 
mun. 20, 425 (1976); T. Ando, Z. Phys. B 26, 263 (1977); 
Advances in Solid State Physics, vol. 48, edited by R. 
Haug (Springer, Berlin 2009). 

[7] For a recent review see B. Spivak, S.V. Kravchenko, S.A. 
Kivelson, and X.P.A. Gao, arXiv:0905.0414 

[8] M. Zaluzny, Phys. Rev. B 47, 3995 (1993); K. Craig, 
et al. Phys. Rev. Lett. 76, 2382 (1996); S. Lutgen, et 
al, Phys. Rev. B 54, R17343 (1996); I. Shtrichman et 
al. Phys. Rev. B 65, 035310 (2001); T. Miiller, W. Parz, 

G. Strasser, and K. Unterrainer, Phys. Rev. B 70, 155324 
(2004). 

[9] M. I. Stockman, L. N. Pandey, L. S. Muratov, and 
Th. F. George, Phys. Rev. B 48, 10966 (1993); 
A. A. Batista et al, Phys. Rev. B 66, 195325 (2002); 

H. Wijewardane and C. A. Ullrich, Appl. Phys. Lett. 84, 
3984 (2004); J. H. Li, Phys. Rev. B 75, 155329 (2007). 

[10] E. Collin et al, Phys. Rev. Lett. 89, 245301 (2002). 
[11] P. M. Platzman and M. 1. Dykman, Science 284, 1967 

(1999) ; S. A. Lyon, Phys. Rev. A 74, 052338 (2006). 
[12] D. Konstantinov et al, Phys. Rev. Lett. 98, 235302 

(2007); D. Konstantinov, Yu. Monarkha, and K. Kono, 
Low Temp. Phys. 34, 377 (2008). 

[13] R. S. Crandall and R. Williams, Phys. Rev. A 5, 2183 
(1972); M. W. Cole, Rev. Mod. Phys. 46, 451 (1974); 
Yu. P. Monarkha and V. B. Shikin, Sov. J. Low Temp. 
Phys. 8, 279 (1982). 

[14] M. Saitoh and T. Aoki, J. Phys. Soc. Jpn. 44, 71 (1978). 

[15] D. Ryvkine, M. J. Lea, and M. L Dykman, 
Abstracts of the APS March Meeting (2006), 
http://meetings.aps.org/Meeting/MAR06/Event/45165i 
Invited Presentations at the Meeting on "Floating Elec- 
trons on Helium for Quantum Computing", Paris (2006), 
http:/ /www.princeton.edu/^lyon/paris_2006 

[16] P. Glasson et al, Physica E 22, 761 (2004). 

[17] C. Fang- Yen, M. L Dykman, and M. J. Lea, Phys. Rev. 
B 55, 16272 (1997). 



5 



SUPPLEMENTARY MATERIAL: the talk by M. I. Dykman at the International Meet- 
ing on "Floating Electrons on Helium for Quantum Computing", Paris (2006) 



HYSTERESIS AND BLEACHING OF 
ABSORPTION BY ELECTRONS ON HELIUM 

D. Ryvkine,^ M.J. Lea,^ and M.I. Dykman^ 

^ Department of Physics and Astronomy, Michigan State University 
^ Royai Hoiioway, University of London 



•Dynamics for slow energy relaxation 
•Absorption bleaching 
•Many-electron hysteresis 
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Microwave absorption and saturation 



Stark-shift transition frequency by a field to tune to 1- 2 resonance 
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Interpretation: absorption saturation in a 
two-level system 



Conventional absorption saturation 
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Nearly resonant driving: 

//^ =-F2^J 2)^^(1 1 exp(-/^^0 + h.c. F = \eE^^z 
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Weak field absorption: (T = ;[flFI^, Z = ^U^o + (^21 



Weak-field linewidth is = F + P. 



Power absorbed per unit time is aco^ 



2a/r 



Strong field absorption: a = rz \ F \^ /(T + 2z \ F \^) 



Electrons: bands of in-plane motion instead of energy levels 
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Electron relaxation mechanisms 



Electron-electron coupling: thermalization of in-plane motion over time CD , 



Electron in-plane momentum distribution p^^ (p) oc exp(— / 2mT^ ) {k^ = 1) 




One-ripplon/vapor atoms scattering is quasi-elastic 

1/2 

and short-wavelength, q » Ifl 

' intraband scattering, T"" 

• interband scattering, F^^ ((2 ^ P) 
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Energy relaxation: two-rippion/phonon scattering 



Master equation for many-electron system 



Field Hamiltonian for spatially uniform resonant radiation F = ^eE^^Zi2 

=-^Zj2L(l|expH6>fO + h.c. 
Frequency detuning is small: Sco = co ^ - (s^ - s^), \Sco\« cd^ 

2 3/2 1/2 

Short-wavelength scattering, > co^ = {27ie n I m) 

effectively single-electron kinetic equation for a strongly correlated electron system 

In the band index representation 

X 

^ T. Ando, 1978 



Detailed balance for fast in-plane thermalization T^^ = T^^ exp[(6*^ ~^B)^^e^ 



Energy relaxation 



Energy balance equation 



dt 



l\/licroscopic mechanisms: 

• Energy diffusion from one-ripplon scattering 

• Two-ripplon scattering, I + q2 I « I q^ 2 I 

• Decay into phonons: modulation of the He dielectric 
constant, Qu « 
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Two-ripplon scattering, kinematic 
coupling: 



Distribution over bands 



Slow energy relaxation, F T » 1 : the rate of field induced transitions does not 
have to beat the rate of 2^1 relaxation. Excited states are populated and absorption is 
bleached already for x\^\ « due to electron heating, x — ^0 1(^1 ^ ^^^) 

Reminder: absorption saturation for a two-level 
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Equation for electron temperature 
energy is counted off 6^ 
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x\FU<r-' 



thermal distribution 
over bands in the 
stationary regime 
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Absorption decrease I 




T —T I cOj, 

e e F 



Absorption decrease II 
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Many-electron line shift 



Strongly correlated electron liquid for (;r /l)^^^ » 

Different distance from He surface in different states 
leads to dependence of electron transition frequency 
on states of neighboring electrons: 
many-electron Stark shift 

"Many" nearest neighbors: mean-field approximation 
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Energy balance equation 
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Many-electron bistability 



Self-consistent equation for ee temperature 

T y CT ^ T 
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Z(f^) has a narrow peak for Set) = fl(T*) if 
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Constant energy relaxation rate approximation 




weak field 



strong field bistability 
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Many-electron bistability 



Self-consistent equation for ee temperature 

I^^ = \F[x{f,)T, f^=^ 



\-e 



Z(fj has a narrow peak for Sco = ^^(7^/) if 
[dQ/dfj 

Kinematic two-ripplon scattering 
weak field 
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Strong field bistability 
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Conclusions 



> Absorption saturation is accompanied by bleaching from 
electron heating. 

> Many-electron shift of transition frequency leads to absorption 
hysteresis for low helium temperatures 

> Electron energy relaxation can be studied via nonlinear 
absorption experiments 



